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ABSTRACT 



We report here the first evidence for planetesimal infall onto the very young Herbig Be star 
LkHa234. These results are based on observations acquired over 31 days using spectroscopy 
of the sodium D lines, the He I 5876A, and hydrogen Hq lines. We find Rcdshiftcd Absorp- 
tion Components (RAC) with velocities up to 200 km/s and very mild Blueshifted Absorption 
Components (BEC) up to 100 km/s in the Na I lines. No correlation is observed between the 
appearance of the Na I RAC & BEC and the H^ and He I line variability, which suggests that 
these (Na I RAC & BEC) are formed in a process unrelated to the circumstellar gas accretion. 
We interpret the Na I RAC as evidence for an infalling evaporating body, greater than 100 km 
in diameter, which is able to survive at distances between 2.0 to 0.1 AU from the star. The 
dramatic appearance of the sodium RAC and mild BEC is readily explained by the dynamics of 
this infalling body making LkHa234 the youngest (age ~ 0.1 Myr) system with evidence for the 
presence of planetesimals. 



Subject headings: stars; formation 
vidual (LkHc234) 

1. Introduction 



stars: pre-main-sequence — circumstellar matter — stars: indi- 



LkHa234 (mi,=11.9mag) is a star embedded in 
a nebula associated with the NGC 7129 star form- 
ing region and is at a distance of about 1250 par- 
sees. This is one of the youngest known Herbig 
Be objects and has an age ~ 0.1 million years 
and spectral type of B5 (Fuente et al. 2001). 
Herbig Ae/Be stars are intermediate mass (2 to 
IOM0) pre-main-sequence stars (Herbig 1960). A 
number of Herbig Ae stars have been found to 
possess circumstellar disks and show both photo- 
metric and spectroscopic variability (Herbst et al. 
1999). LkHQ,234 though much hotter and more 
massive also exhibits these properties (Fuente et 
al. 2001, Herbst et al. 1999, Polomski et al. 2002). 
The star shows photometric variability of 1.2 mag- 



^ Based on observations from the 9.2m Hobby-Eberly 
Telescope 



nitude in the V band and is known posses a cir- 
cumstellar disk with the disk mass to the stellar 
mass ratio (M^/M*) < 0.02. A detailed descrip- 
tion of LkHQ234 and its circumstellar region can 
be found in Fuente et al. (2001, and references 
therein). 

Spectroscopic variations in metallic lines (like 
Na I, Ca II, etc.) in young A type stars like fi 
Pictoris (20 to 100 million years) are often inter- 
preted as the signature of star grazing comets or 
falling evaporating bodies with sizes of tens of kilo- 
meters (Grinin et al. 1996, Smith & Terrile 1984, 
Thebault & Beust 2001). The only Herbig Be star 
where grazing cometary transient events have been 
inferred so far is HD 100546 (Viera et al. 1999). 
However, HD100546 is much older (>10 Myrs) and 
cooler (spectral type ~B9) compared to LkHa234. 
High resolution spectroscopic monitoring of very 
young (0.1 to 1 Myrs) Herbig Be stars are scanty 



and uncharted. 

Here we report dramatic variations in the spec- 
troscopic lines of Na I 5890A(D2), and 5896A(D1), 
He I 5876A, and of LkH„234 over a period of 
31 days during Oct.-Nov. 2003. The observations 
consists of high resolution spectra (R=30,000) 
from the Hobby-Eberly Telescope (HET) using the 
High Resolution Spectrograph (HRS). Section 2 
describes observations and in section 3 we discuss 
the results and section 4 gives the conclusion. 

2. Observations 

We have obtained high signal to noise ratio 
and high resolution spectra (R=30,000, '~0.2A) of 
LkHa234 covering He I 5876A, Na I D2 & Dl lines 
and the Balmer H^ at 6563A. The observations 
were taken with the 9.2m Hobby-Eberly Telescope 
(HET) (Ramsey et al. 1998) using its High Resolu- 
tion Spectrograph (HRS) (TuU 1998). The echelle 
spectral orders were optimized to cover the wave- 
length range between 5700 to 6750A. This allowed 
the spectral lines of interest to be recorded by the 
"blue" CCD of the spectrograph that has fewer 
number of bad pixels than the "red" CCD. The re- 
sults comprise of five data sets over a time period 
of 31 days. The observation dates were: 7th, 13th, 
23rd, 27th October and 8th November 2003. On 
each night of observation, three ten minutes ex- 
posures on LkH„234 were taken which were com- 
bined to get high signal to noise ratio (S/N = 70 to 
100). All data were reduced using standard IRAF 
routines under ECHELLE tasks. Every spectrum 
was bias subtracted, and flat-fielded. Spectral Cal- 
ibrations were done using Thorium-Argon lamp 
spectra taken either immediately before or after 
the observations. We found that the wavelength 
calibration is accurate up to O.OSA. 

H„ (6563A), He I 5876A, and Na I D2 & Dl 
spectral lines arc shown in figures 1 and 2 re- 
spectively and the dates of observations are noted 
along side of the spectra. The x-axis is in terms of 
relative velocity (km/s) with respect to the center 
of the lines in the stellar rest frame and the y-axis 
shows line intensities normalized with respect to 
the stellar continuum. In figure 2 the zero of the 
relative velocity scale is with respect to the stellar 
rest frame of the sodium D2 line. 



3. Results &: Discussion 
3.1. Gaseous accretion 

We illustrate in figure 1 the variability in Hq, 
and He I lines. These are seen as stellar photo- 
spheric absorption lines in the spectra observed 
on 7th Oct. 2003. However, in all the later days 
of observations the H^ appeared to be in emis- 
sion consisting of two components, the Red shifted 
Emission Component (REC), and the Blue shifted 
Emission Component (BEC). The Hq, equivalent 
width varied from +3A(on 7th Oct.) to -50A(13th 
Oct.) in just 6 days and it continued to rise up to 
-100A(27th Oct.) before it decayed to -74A(8th 
Nov.). The relative velocity between the REC and 
BEC varied between 150 km/s to 225 km/s with 
the REC stronger than the BEC by a factor of 2.5 
to 3.5. 

The spectra obtained on 7th Oct. 2003 do not 

show any sign of emission or activity in Hq, or 
He I. This observation shows the star in a quies- 
cent state, with neither accretion nor stellar wind 
affecting the absorption line profile. The varia- 
tions in the He I line are found to correlate with 
the variations in the Hq,. Whenever the H^ BEC 
and REC are stronger the He I either fills up due 
to mild emission or shows a very large absorp- 
tion width (almost 1.7 times the photospheric line 
width). The line width varied from 150 km/s to 
260 km/s. A B5 star does not have sufficiently en- 
ergetic photons to ionize He I (Jaschek & Jaschek 
1995), so the filling up of the He I absorption line 
due to mild emission is an indicator of accretion 
phenomena (de Winter et al. 1999). We find the 
disk luminosity to be about 1% of the stellar lu- 
minosity (L« ^IO^Lq) with maximum disk tem- 
perature of ~21000K, if we consider an accretion 
rate of lO-'^M© per year (Natta et al. 2000). 

Therefore, the double peak H„ (REC & BEC) 
profiles and the He I mild emission are because of 
orbiting accreting gas on to the star and as well 
as outfiow of circumstellar gas (BEC in Hq due to 
accretion heat). However, when the accretion lu- 
minosity is not strong enough to produce enough 
ionizing photons, the He 1 line appears as a very 
broad absorption feature. The broad absorption 
feature is the effect of the combination of the stel- 
lar photospheric line and various unresolved cir- 
cumstellar gas components corresponding to REC 
and BEC of Hq. 



2 



The effect of gaseous accretion are probably 
also seen on the nebular emission lines. We found 
a correlation between the total equivalent width 
in Ha and the peak nebular emission in Hq and 
Na I. Figures 1 & 2 show that the nebular Hq, 
and Na I lines were of minimal intensity (~1.06 
above the continuum) on 7th Oct. The same lines 
are strongest (£[„ ~7.9, on the blue wing of the 
REG, and Na I D2 line ~3.0) on 27th Oct. when 
the Hq, equivalent width is maximum (-lOOA). We 
plan to model the impact of the episodic accretion 
on the nebula using a photo-ionizing code (for eg. 
the Cloudy code) which will be published along 
with the entire spectra of LkHQ234 elsewhere in 
the near future. 

3.2. Planetesimal infall 

Figure 2 shows the variability in the Na I D2 
and Dl lines. Redshifted Absorption Compo- 
nents (RAC) were observed in Na I lines on only 
one night (13 Oct. 2003) with a maximum red- 
shifted velocity of 200 km/s. The RAC's associ- 
ated with both the Na I lines arc of similar depth 
(for redshifts < lOOkm/s) which is an indicator 
of saturation (unshielded Na I column density = 
10^^/cm^)(de Winter ct al. 1999). However, the 
maximum depth of these components is about 30% 
of the stellar continuum implying partial coverage 
of the stellar photosphere. In addition to Na I 
RAC, very mild BEC (60 to 100 km/s and 1.05 
to 1.07 above the continuum) are also seen on the 
13th Oct. spectra. The variations in He I and 
Ha are observed even when no Na I RAC & BEC 
are seen (see figures 1 & 2). Thus the appear- 
ance and disappearance of the Na I RAC & BEC 
are uncorrelated with the variations in H^ and He 
I line profiles. We conclude that we witnessed a 
transient phenomenon on 13th Oct. 2003 whose 
effects lasted a few days at most. 

Prom the Keplerian dynamics of an infalling 
object (Beust et al. 2001) the most redshifted 
Na I absorption component velocity (200 km/s) 
should correspond to a distance of 0.1 AU from 
the star. At this distance both the number of 
ionizing photons and intensity of the stellar wind 
will be too high for any unshielded Na I to sur- 
vive (SorcUi ct al. 1996). Their calculations using 
photo-ionizing codes under local thermodynamical 
equilibrium (LTE) and non-LTE show that only a 
solid body like a comet or asteroid can approach a 



hot star this close before it disintegrates and evap- 
orates. Further, theoretical model calculations for 
stars up to B9 show that the size of solid bod- 
ies which arc able to survive up to a distance of 
0.1 AU should be at least 100 km in diameter 
(Beust et al. 2001). The mild Na I BEC could 
then be part of the falling evaporating body be- 
ing blown away by the stellar wind or radiation 
pressure. However, this will also mean that the 
falling solid body may have a high eccentric orbit 
so that the blown away material is not projected 
against the stellar surface. Our Na I RAC and 
BEC observations arc consistent with such an ob- 
ject dissociating between 2 AU and 0.1 AU from 
the B5 star LkH„234. 

In many Herbig Ae stars, variations in Na I 
RACs and BAC's usually correlates with H^ and 
He I variations (de Winter et al. 1999, Grinin et al. 
1996) . Such Na I variations are now thought to be 
generated by magneto-hydrodynamic funnelling of 
neutral gas on to the star rather than solid body 
infall (Beust ct al. 2001, Mora et al. 2002, Natta 
et al. 2000). No such correlation is present in the 
current LkHa234 data set. Of particular signifi- 
cance is the observation on 27th Oct. 2003 when 
the REG and BEC of H„ are the strongest. The 
absence of any Na I RAC and BEC in this observa- 
tion is a strong argument against the Na I line vari- 
ations being caused by neutral matter funnelled 
onto the star by magnetic fields. The variability of 
Ha and He I in LkHc234 are due to episodic gas ac- 
cretion, while the observed variations in Na I lines 
is a transient event that exhibits the dynamical 
signature of a body (>100 km in diameter) falling 
onto the star. It is worth mentioning here that 
recent studies on meteors from cometary origin by 
Trigo-Rodriguez et al. (2004) have shown greater 
sodium abundances than those expected for inter- 
planetary dust particles and chondritic meteorites 
and Potter et al. (2002) have detected comet like 
sodium tail from planet Mercury, and also the dis- 
covery of presence of sodium tail in comets (Cre- 
moncse et al. 1997). Thus an infalling solid body 
of asteroid size and deficient in H can produce the 
observed Na I RAC and BEC and no correlation 
with the H and He I lines. 

The dusty circumstellar disk of the young star 
LkH„234 could be a newly formed protoplanetary 
disk consisting of solid bodies of different sizes 
(Fuente et al. 2001). Theoretical models of for- 
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mation of planetary bodies in circumstellar disks 
around a IM© mass star show that planetesimals 
of sizes of up to a few hundred kilometers can 
quickly form within the first 10"'' years and the typ- 
ical cumulative number of such bodies can be 10^ 
to 10^ between 0.5 to 1.5 AU from the star (Lis- 
sauer 1993, WeidenschiUing ct al. 1997). While 
such numbers are not known for B type stars, 
LkH„234 (age ~ 0.1 Myrs) may have high fre- 
quency of infall events (Grady et al. 2000). 

4. Conclusion 

We have presented here high resolution spec- 
troscopic monitoring of a Herbig Be star LkHQ,234 
in the Na I (D2&D1), He I (5876A) and H„ lines. 
We found no correlations between the Na I and 
the He I and Hq, line variability in the five data 
sets over a period of 31 days (7th Oct. 2003 to 
8th Nov. 2003) though we do find a correlation 
between He I and H^. Thus the origin of vari- 
ation in the Na I line (RAG up to 200 km/s and 
BEG up to 100 km/s seen only on 13th Oct. 2003) 
is different from those of H^ and He I. While the 
Hq, and He I line variations are due to episodic 
gaseous accretion, the Na I RAG indicate a dra- 
matic transient event of solid body infall on to the 
star and the BEG as the blown away parts of the 
solid body not projected against the stellar sur- 
face. Gonsidering the Keplerian dynamics and the 
harsh environment of a B5 star, we estimate that 
a solid body of size >100 km broke up and disin- 
tegrated at a distance between 0.1 to 2.0 AU from 
the star. This makes LkHQ,234 the youngest sys- 
tem (~0.1 Myrs) with evidence for protoplanetary 
bodies of asteroidal size. We plan to pursue fur- 
ther spectroscopic monitoring of the star with the 
HEX to determine the frequency of such events 
and further understand the complicated dynamics 
of this very young circumstellar environment. 
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Fig. 1. — The left panel shows the spectra of 
LkHa234 around the Ha(6563A) emission line and 
the right panel around Hel (5876 A). The x-axis is 
in relative velocity with respect to the rest wave- 
lengths of Hq and Hel respectively. The y-axis 
is the flux normalised with respect to the stellar 
continuum. The dates of observations are shown 
with each spectrum. The narrow 20 km/s width 
nebular emission line is seen near the; core of 
the stellar photospheric absorption line on 10-07- 
03 and so are many nebular absorption lines. The 
nebular Hq line can also bo seen on the blue wing 
of the REC on the other dates of observations. 
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Fig. 2. — The spectra of LkHa234 in the region of 
Na I D2(5890A) and D 1(5896 A) absorption fea- 
tures. The x-axis has been converted into relative 
velocity (km/s) with respect to the Na I D2 line 
at the stellar rest frame. The y-axis is same as 
in figure 1. The Red shifted Absorption Compo- 
nents (RAC), the nebular emission and the inter- 
stellar absorption linos are marked. Tlio dates of 
observations are shown along with each spectrum. 
Also seen are the very mild Blueshifted Absorp- 
tion Components (BEC; on 10-13-03) on the blue 
side of the Na I (D2&D1) interstellar absorption 
lines. The interstellar Na I lines are offset by 27 
km/s from the rest frame of the star. 
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